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Abstract: From various cosmological, astrophysical and terrestrial requirements, we de- 
rive conservative upper bounds on the present-day fraction of the mass of the Galactic 
dark matter (DM) halo in charged massive particles (CHAMPs). If dark matter particles 
are neutral but decay lately into CHAMPs, the lack of detection of heavy hydrogen in sea 
water and the vertical pressure equilibrium in the Galactic disc turn out to put the most 
stringent bounds. Adopting very conservative assumptions about the recoiling velocity of 
CHAMPs in the decay and on the decay energy deposited in baryonic gas, we find that 
the lifetime for decaying neutral DM must be > (0.9 — 3.4) x 10 3 Gyr. Even assuming the 
gyroradii of CHAMPs in the Galactic magnetic field are too small for halo CHAMPs to 
reach Earth, the present-day fraction of the mass of the Galactic halo in CHAMPs should 
be < (0.4 — 1.4) x 10~ 2 . We show that redistributing the DM through the coupling be- 
tween CHAMPs and the ubiquitous magnetic fields cannot be a solution to the cuspy halo 
problem in dwarf galaxies. 



Keywords: dark matter - galaxies: haloes - galaxies: kinematics and dynamics 



galaxies: magnetic field; 



*E-mail:jsanchez@astroscu.unam.mx 



Contents 



|]. Introduction |] 

{2L Assumptions |2] 

j3j_ Pregalactic constraints |3| 

|[ The density of CHAMPs in the Galaxy | 

|4. 1| Shielding the disc with magnetic fields || 

4.2 Energy gain and loss of CHAMPs in the disc ^ 

[4.3| Constraints from sea water searches || 

[5]. The global magnetic support of halo CHAMPs 10 

§. Ram pressure stripping and collisions of galaxy clusters [l4j 

0. Conclusions 15 

|A] . Recombination and ionization of CHAMPs in the Galactic corona 16 

|B|. Solving the differential equation for n\ d 17 



1. Introduction 

So far, other than its gravitational interaction, the detailed properties of the dark matter 
(DM) are still largely unknown. The lack of direct detection suggests that DM particles are 
stable, neutral and weakly interacting. However, it is important to know how well these 
properties are constrained from an observational point of view. In recent years, mainly mo- 
tivated by the cuspy problem of dark haloes in low-surface brightness galaxies and dwarf 
galaxies, and the excessive abundance of satellite galaxies inferred in cosmological simula- 
tions, much work has been done to study the cosmological and astrophysical implications 
of other variants, such as decaying, collisional or annihilating DM. New interesting phe- 
nomena at galactic scales, ignored under the assumptions of collisionless and neutral DM, 
arise if a fraction of DM is made up by massive particles with electric charge (CHAMPs). 
Several theoretical physics models beyond the Standard Model (SM) have shown the possi- 
bility of CHAMPs (e.g., @, §). The existence of CHAMPs is well-motivated in the model 
of super Weakly Interacting Massive Particle (super- WIMP) dark matter Q. In super- 
symmetric theories, some popular examples of CHAMPs include a slepton, such as the 
super symmetric staus (e.g., [||, |]]). 
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The hypothesis that the DM is made up by a mixture of bare CHAMPs and neutra- 
CHAMPs (a neutral bound atom formed by a CHAMP of electric charge —1 and a proton) 
was considered in the late eighties |J. Dimopoulos et al.|| noticed that the astrophysical 
and terrestrial limits are hardly compatible with such a scenario (see j7j and [|| for updated 
reviews). Nevertheless, all these constraints were derived for the standard flux of particles 
at Earth from the Galactic halo, assuming it was constant over time. If DM particles are 
originally neutral and decay lately into charged particles, many constraints can be avoided 
and one must reevaluate earlier bounds. 

This paper is organized as follows. In section ^j, we describe our assumptions, which 
were especially designed to minimize the potentially disastrous effects of charged DM. 
The constraints from big bang nucleosynthesis (BBN) and cosmic microwave background 
(CMB) anisotropy on the lifetime of decaying neutral DM particles are discussed in section 
||[ In sections |] and [|, we examine the distribution of CHAMPs in the disc and in the halo 
of the Galaxy. We will discuss the implications of the lack of detection of anomalously 
heavy hydrogen in sea water and the vertical (magneto-)hydrostatic configuration of the 
Galactic disc on the fraction of halo CHAMPs. Other physical implications of CHAMPs 
embedded in galaxy clusters are briefly discussed in section ||. Concluding remarks are 
given in section [j]. 



2. Assumptions 

In order to remedy purported problems with the collisionless CDM family of cosmological 
models on galactic and galaxy cluster scales, Chuzhoy & Kolb [||] have revived the possibil- 
ity that a significant fraction of the DM in haloes is made up by unneutralized CHAMPs. 
These authors claim that the distribution of CHAMPs in galaxies may be altered by the 
coupling between CHAMPs and ubiquitous magnetic fields. For instance, CHAMPs may 
be depleted from the central parts of galaxies, erasing the DM cusp, if they are accelerated 
through the Fermi mechanism in supernova shocks. 

Here we explore a generic scenario in which neutral dark matter, denoted by x, decay 
with lifetime Td ec into non-relativistic and exotic massive particles with electric charge. The 
model of Chuzhoy Sz Kolb |J corresponds to a scenario where the decay lifetime is shorter 
than the age of the Universe. The decay of DM into another dark or SM species could have 
a bearing upon possible problems with the ACDM scenario, e.g., the reionization of the 
Universe, the structure formation at small-scales and the low abundance of satellite galax- 
ies, the formation of galactic cores, the synthesis of light elements in the early Universe, 
the origin of ultra-high energy cosmic rays, the positron excess observed by PAMELA, the 
Tully-Fisher relation with z or the gas fraction of galaxy clusters. However, we will focus 
on the exclusion range of parameters. 

Stringent limits on DM decay into photons or SM particles have been derived from 



diffuse 7-ray observations |10[], the effects on BBN [11] and from the reionization history of 
the Universe (e.g., fl2| , fOJ). In order to set an upper limit on the fraction of CHAMPs as 
generous as possible, we will consider the most favourable and simplest scenario to permit 
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the maximum amount of charged particles in galactic haloes. Our model starts from the 
following, rather artificial, assumptions: 

(1) The coupling strength of the decay of a neutral DM particle into two (electrically) 
charged particles, X± and X^ , plus a very light (or massless) weak interacting particle X v , 
is large. Therefore, x dominantly decays without the emission of photons or Z bosons. The 
X-decay mode is thought to be analogous to the neutron /3-decay but in the dark sector. 
Note that is not the anti-particle of X^ and thus they may have different masses. To 
keep the discussion manageable, however, we will assume that both particles have the same 
mass. 

(2) X^~ and X^ are stable and cannot decay. This can be accomplished if they are the 
lightest particles carrying conserved 'dark baryon' number B' and 'dark lepton' number 
L', respectively. Conservation of these quantum numbers in the decay is met if % and X^~ 
are dark baryons and X^ is a dark lepton and X u a dark anti-lepton. 

(3) The decay \ ~^ ^1 + ^2 + X v is the dominant way through which CHAMPs can 
be produced. This non-thermal production of particles is present in some schemes such 
as SUSY (e.g., 0]). In practice, it is assumed that DM is neutral before freeze-out and 
afterwards it decays into CHAMPs. 

(4) The charged X particles from the decay are non-relativistic so that they are cold 
DM by the time of structure formation and its contribution to the reionization of the 
intergalactic medium does not contradict CMB data. This is fulfilled in models where 
Am/m x <2x 10 -4 , where Am is the mass difference between the initial and final states. 
Although this fine-tuning is unlikely, it is not so rare in nature. For instance, the corre- 
sponding ratio in the classical neutron /3-decay is comparable (~8x 10 -4 ). 

We will see that, even under these somewhat artificial assumptions, the fraction of dark 
matter that can be made by CHAMPS is vanishingly small for astrophysical phenomena 
to be affected and the standard Cold DM cosmology is recovered. The proposed scenario 
is in some ways reminiscent of models already discussed in the literature but introduce 
important conceptual differences. In previous works, decaying DM was introduced to 
dissolve the central cusp in DM haloes and the overabundance of satellite galaxies by the 
depletion ]l5|, |l6[ or energy release in the decays [ 17 , . As suggested by Chuzhoy & Kolb 
§, in the CHAMP model, even if the recoiling velocities of CHAMPs were very small, they 
may be ejected from the central parts of the galaxies by Fermi acceleration in shock waves, 
or from the Galactic disc, making them very evasive for direct terrestrial detection (see Q 
for a discussion). In the next section, we discuss the limits on Td cc imposed by BBN and 
CMB. 



3. Pregalactic constraints 

The standard BBN theory has been well established to predict precisely the primordial 
light element abundances and constrain the number density of long-lived CHAMPs at 
t < 10 5 s. Since the recombination of CHAMPs with protons, a-particles, electrons, or 
other CHAMPs to form neutral atoms, occurs well after BBN, CHAMPs remain bare at 
BBN. If X2 particles, with masses below the weak scale, are present at the BBN and they 
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do not decay into other particles, excessive production of 6 Li and 7 Li may occur only if the 
fractional contribution of negative CHAMPs to the present critical density, Qx, is larger 
than a certain maximum value f2^ ax = 3 x 10 -6 (e.g., [19, pC] ], and references therein) 1 . 
For masses above the weak scale, this bound can be weakened. In our case, the abundance 
of CHAMPs increases in time due to % decays and this constraint only applies at t < 10 5 
s. By imposing that, at most, a particle number fraction of f2™ ax /f2 c of x's, with Q c the 
present density of cold DM (£l c = 0.23), has decayed by the end of BBN, t ~ 10 5 s, we 
obtain r^cc > 200 yr. 

A much more stringent constraint can be placed by studying the effect on the CMB 
anisotropy of a scattering interaction between charged CHAMPs and the photon-baryon 
fluid. Throughout the paper, we will refer to "charged CHAMPs" to indicate both the 
free CHAMPs and the bound states of CHAMPs having a net electric charge. Kohri &: 
Takahashi [21] have shown that most of the negative CHAMPs are captured by a particles, 
forming charged CHAMPs. Hence, all the CHAMPs are expected to be in a charged state at 
the epoch of recombination. Since charged particles have low velocities, they can thermalize 
very quickly by their interaction with the baryons [|| . The effect on the CMB of coupling 
baryons with charged CHAMPs is equivalent to a standard model with a larger value of 
Ofo. Given that the WMAP uncertainty in is < 3%, the fraction of DM particles that 
are allowed to decay before recombination epoch, tji, is < 0.03(r2b/f2 c ) ~ 6 x 10 -3 , where 
Ob/f^c is the fraction of baryonic energy density f^/i 2 , relative to that of the (cold) DM, 
S7 c /i 2 . Consequently, CMB power spectra is obtained if Td ec > £r x 10 3 /6 ~ 7 x 10 7 yr, 
where we have used tR ~ 376,000 yr (see ]23] for the case of millicharged CHAMPs). 

For invisible decay to weakly interacting particles such as neutrinos or X u , Gong & 
Chen [23] constrain the decay lifetime to t^ cc > 0.7 x 10 3 £ Gyr, where £ is the fraction of 
the rest mass which gets converted to neutrinos or X v . Our assumption (4) implies that 
£ ~ 2x 10 -4 and, consequently, we should constrain ourselves to models with Td ec > 1.5 x 10 s 
yr. 

In contrast to the model of De Rujula et al. ||, the contribution of the pressureless 
X component, and not neutraCHAMPs, is dominant in driving galaxy formation. In this 
sense, our model is closer to the standard collisionless CDM. 



4. The density of CHAMPs in the Galaxy 

After recombination, when photons and baryons become noninteracting, baryons and 
charged CHAMPs fall into the gravitational wells formed by pressureless x-particles and 
neutraCHAMPs. In the halo of a galaxy like ours, the temperature becomes high enough 
to ionize hydrogen and superheavy hydrogen (X^e~) -see the Appendix |A|-. In contrast, 
although the equilibrium fraction of (X^p) in the halo is small, neutraCHAMPs formed 
at pregalactic stages (we will refer to them as 'primordial neutraCHAMPs) can survive a 
long time compared to the Hubble time, without being dissociated (see the Appendix |A|). 



Nevertheless, the initial density of primordial neutraCHAMPs is expected to be small [21 



lr rhe CHAMP-to-entropy ratio Y x was converted to Q. x using Q x h 2 = 2.73 x 10 11 Y x (m x /1 TeV). 
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and they might be converted to (X^ct) ions by charge-exchange scattering and then they 
behave as charged CHAMPs in the lifetime of the Galaxy [||. 

In order to constrain the abundance of CHAMPs in the Galactic disc and in the halo 
we must consider the following processes: (1) magnetic fields can prevent the charged 
CHAMPs in the halo to penetrate the Galactic disc, (2) charged CHAMPs may be ejected 
from the disc and blown either back to the halo, or right out of the galaxy if charged 
CHAMPs are shock accelerated by supernovae (e.g., [||, (3) Neutral \ particles have 
no difficulties to penetrate the disc and may have a decay when they are crossing the disc, 
replenishing the disc with fresh CHAMPs. In the following, we study conditions for which 
the abundance of CHAMPs is small in the disc, making its detection unlikely, but large in 
the halo. 



4.1 Shielding the disc with magnetic fields 

In principle, a charged CHAMP may lose all its kinetic energy by interacting with the 



interstellar gas as it intersects the Galactic disc [24|. However, the penetration of unneu- 
tralized coronal CHAMPs along the Galactic disc is impeded by the presence of Galactic 
magnetic fields. Chuzhoy & Kolb || claim that charged CHAMPs cannot cross the disc if 
they have masses mx < 10 8 TeV. A more precise calculation taking into account that the 
Galactic magnetic field is not plane-parallel is given below. 

It is well-known that when charged particles interact with a magnetized body, a bound- 
ary layer that divides two regions with different conditions is created pEf |. Thus, charged 
particles will penetrate this boundary by some distance before they are turned around by 
the v x B force (Figure ||). 
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Boundary layer 



Figure 1: Sketch of the structure of the magnetic boundary layer formed by the partial penetration 
of the charged particles before they are deflected back. We also illustrate the corresponding orbits 
of particles in the neighborhood of such boundary for both types of charge, that is, negative and 
positive. The degree of shielding depends on the topology of the magnetic field. In galactic discs, 
the magnetic field lines are not plane-parallel. 

The boundary layer is formed because of the partial penetration of the charged particles 
before they are deflected back. A schematic representation of the orbits described by 
negative and positive charged particles in the neighborhood of the magnetic boundary are 
drawn in Fig. [l|. 

The degree of shielding by magnetic fields depends on the configuration of the magnetic 
field. In the case of the magnetized Galactic disc, the magnetic field is not plane-parallel 
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and the propagation of charged particles is more complex because the magnetic field has 
a turbulent component, i.e. B = Bq + b, where Bq is the regular (homogeneous) magnetic 
field and b denotes the turbulent field. The particle motion is determined not only by the 
average magnetic field but also by scattering at field fluctuations, a stochastic process which 
requires the solution of transport equations with particle ensembles. Particle propagation 
in turbulent fields can be understood as a diffusive process, reason why we consider the 
spatial diffusion of halo CHAMPs into the galactic disc. 

As occurs when one considers the escape timescale of cosmic rays from the Galactic disc 
(e.g., |2q| ), the timescale for the penetration of halo CHAMPs is governed by the diffusion 
timescale across the galaxy disc thickness, r^iff. The ordered magnetic lines follow the spiral 
pattern, come out of the Galaxy disc and unfold in the halo. However, due to the presence 
of the tangled (turbulent) component of the magnetic field, the penetration timescale of 
halo particles through the magnetic spiral arms is much longer than r^m because they must 
diffuse a distance much longer than the galaxy disc thickness H. The diffusion timescale 
Tdiff across H for a halo CHAMP, is bracketted in the range: 

tt2 tt2 

' <r diS <—, (4.1) 



2D|, ^ ' mtt ~ 2ZV 

where Dy and D± are the diffusion coefficients parallel and transverse to the mean com- 
ponent of the magnetic field, which is observed to be parallel to the disc. 

The magnetic field that feels a charged CHAMP moving within the disc can be con- 
sidered static because Alfven waves propagate with velocities of the order of the Alfven 
speed v a ~ 6 km s , which is smaller than the typical velocities of CHAMPs > V3cr v , 
where a v ~ 150 km s _1 is the one-dimensional velocity dispersion for halo particles. The 
diffusion coefficients depend on the turbulence level r] = (1 + (-Bo) / (^ 2 )) _1 > anc ^ on the 
rigidity x = 27rrL/A max , with rx the Larmor radius defined with respect to the total mag- 
netic field and A max the maximum scale of the turbulence ~ H/2 [26, Observations 
of the Galactic polarized synchrotron background yield 1 < (6 2 ) /(-Bg) < 9 (^|, and 
references therein), implying that 0.5 < r/ < 0.9. Since t^is scales as the inverse of the 
diffusion coefficients and those are essentially a monotonic function of r/, we use r\ ~ 0.5 in 
our estimate of D± in order to give an upper limit on the diffusion timescale. From the 
numerical result by [p6[ , we know that D±/(rLv) ~ 0.3 for Kolmogorov turbulence with 
r] = 0.5 and x between 0.05 and 0.4, we find that the diffusion timescale average over the 
velocity distribution is 

/ / 5 # 2 \ ^ ( H \ 2 f m x \-i 
Tdiff -\3^/ = 1 ° Gyr l3^cJ KWfeV) 

( — .W-g-V (4.2) 



V 150km s -1 / \5nG 

The values of the turbulence level and rigidity depend on the galactocentric distance but 
also on the azimuthal angle in the disc, because both the regular and the turbulent fields are 
commonly more intense within the spiral arms ( , and references therein) . The efficiency 
of the magnetic shielding can be reduced along Galactic magnetic chimneys. Moreover, it 
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is likely that CHAMPs are accelerated to much higher velocities by supernova shocks as 
soon as they penetrate inside the disc, decreasing r^g further. In the most optimistic 
scenario where all these effects can be ignored, the present configuration and strength 
of the Galactic magnetic field can prevent diffusion of (unaccelerated) charged CHAMPs 
across the Galactic disc in the lifetime of the disc for mass particles rax < 10 6 TeV. Note 
that the corresponding gyroradius for a mass of 10 6 TeV moving at 300 km s _1 in a field 
of 5/uG is 0.2 pc. The equation that governs the number of CHAMPs in the disc will be 



discussed in section 4.3 



4.2 Energy gain and loss of CHAMPs in the disc 

In the foregoing section we have seen that halo CHAMPs with masses rax < 10 6 TeV may 
have difficulty in penetrating the magnetized Galactic disc, whereas those inside it would 
stay confined to the disc unless they are accelerated. Charged CHAMPs in the disc gain 
energy through electrostatic fields, Fermi acceleration in shock waves, and its descendants 
(e.g., [|(J), and lose kinetic energy due to Coulomb scatterings with electrons and protons 
of the diffuse interstellar gas. 

Consider masses of rax larger than the electron mass 2 . The dissipation timescale due 
to collisions with the electrons is Tdi s = E/\E\, with E = raxv x /2 and 

e 4 

\E\=4:im e — —In A, (4.3) 

ra e v x 

where vx is the velocity of the CHAMP in the interstellar medium, e e is the electron 
charge, and n e is the electron density (~ 0.025 cm~ 3 in the solar vicinity) and the Coulomb 
logarithm has a value of about 20. CHAMPs moving in the Galactic disc may avoid strong 
cooling if the dissipation timescale is greater than the shock acceleration timescale T acc , 



which is > 0.01 Gyr (e.g., [32, ^]). The condition 2rdi s > T acc > 0.01 Gyr implies that 
particles with initial velocities 

vx > v crit ee 150^3 km s- 1 (^fev) ^ ' (44) 

can be accelerated and escape from the disc, whereas those particles with velocities < v CT i t 
are expected to lose kinetic energy until they become neutral by recombining with a proton 
to form a neutraCHAMP (X^p), or an electron to form superheavy hydrogen (X^~e~). 
Basdevant et al. [H] suggested that neutraCHAMPs and superheavy hydrogen in the disc 
will reach thermal equilibrium with the environment and will present turbulent velocities 
as those of interstellar neutral hydrogen, ~ 10 km s . Charge-exchange equilibrium with 
hydrogen dictates that 20% — 40% of superheavy hydrogen in the disc should be ionized. 

Assuming that the velocity distribution of CHAMPs just after the decay of x particles 
is Maxwellian, the fraction F of CHAMPs created in the Galactic disc that will be trapped 
in the disc forming neutraCHAMPs or superheavy hydrogen is 



2 We will not consider the regime mx < rn e because they are excluded for 10 15 < e < 1, where e is the 
electric charge units of e e , the elementary electron charge pifl. 
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After one Hubble time, tjj, the relative abundance of superheavy hydrogen in the solar 
neighbourhood is 



E^l = (F X ID"') (1 - exp [- W ^]) (^^) , (4.6) 

where we have assumed a total density of dark matter of O.OIMq pc~ 3 . Searches for 
interstellar superheavy hydrogen performed by looking at the Lyman j3 absorption in the 
direction of some nearby stars contrain the relative abundance of superheavy hydrogen 



over ordinary hydrogen to < 2 x 10 (e.g., ]24|, 34]). Combining Eqs flO| ) and ( |4.6| ) using 
a v = 150 km s , we find that for Td ec <C tu, the relative abundance is smaller than ~ 10 -8 
provided that mx > 6 x 10 3 TeV. For Td ec — tn, the condition is fulfilled for mx > 4.5 x 10 3 
TeV. 

4.3 Constraints from sea water searches 

If CHAMPs are singly charged, the most stringent bound on the abundance of CHAMPs 
in the disc comes from searches of anomalously heavy sea water 3 . The gyroradius for a 
charged CHAMP of mass 2.5 x 10 3 TeV at 300 km s" 1 is 10 AU in the magnetic field of 
the solar wind (50/uG in the Earth vicinity). Therefore, the arrival of charged CHAMPs at 
Earth cannot be impeded by the magnetic field of the solar wind for mx values in the range 
of interest {mx > 10 3 TeV). The null results of searches of CHAMPs, between 10 3 TeV 
and 10 5 TeV, in ocean water by Verkerk et al. may be used to constrain the admissible 
value for Td ec - 

Denote by n^(R) the number density of positively charged DM particles at galacto- 
centric distance R in the halo. Ignoring the expansion of the dark halo by the recoiling 
velocities in the decay, increases in time by: 

dnf dn y n y 

-rf = --7T = — > (4-7) 
dt dt T dec 

where n x (R,t) is the number density of neutral DM particles in the halo. In the equation 
above, we have neglected the flux of CHAMPs from the disc to the halo since the mass 
of DM in the disc is small as compared to that in the quasi-spherical dark halo. Solving 
this equation, we have n^(R,t) = no(R)[l — exp(— t/T& cc )]. Here, uq(R) is the number of 
X particles at the galactocentric distance R if they would have not decayed. 

As discussed in § |4.2| , for CHAMPs in the disc, we need to differenciate between initially 
slow CHAMPs (v < u crit ) and initially fast CHAMPs (v > v crit ). Slow CHAMPs lose 
energy until they reach thermal equilibrium with the interstellar medium gas, whereas fast 
CHAMPs will be accelerated by supernova shock waves and eventually be ejected from the 
disc. Denote by n^ d the density of neutral and ionized superheavy hydrogen in the 'cold' 
phase (one-dimensional rms velocities of ~ 10 km s _1 ) and by nt d the density of in the 



3 The heating of the interstellar H I gas by the impacts of CHAMPs crossing the gaseous disc imposes a 
constraint of a factor ~ 5 less stringent |B5| . 
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'hot' phase (one-dimensional rms velocities of > 150 km s 1 ). According to our discussion 
in section |4.2| , n^ d increases in time according to 



n+ {t) = Fn (R) 



1 — exp 



t 

T~de< 



(4.8) 



The escape of cold CHAMPs through diffusion is very small and can be neglected. 

For the hot CHAMPs in the disc, we need to consider the evacuation from the disc 
by Fermi acceleration processes, the replenishment of CHAMPs by the decay of neutral 
particles and the flux from the halo to the disc through diffusion. Let r esc the timescale 
for CHAMPs to escape from the disc. The equation for rif, at times t > tj, where is the 
epoch of the formation of the magnetized Galactic disc, is 



dn hd 
dt 



''"dec 



n 



lid 



+ 



n 



n 



lid 



Tdiff 



(4.9) 



with F' = 1 — F. We have made the approximation that the diffusion of CHAMPs from the 
halo to the disc is D±V 2 n + ~ i n h ~ n hd) / TdiS • Using the value of n x previously calculated, 
Equation (fil^) can be immediately solved. For instance, if we define r e fj as 



1 



1 1 

+ : 

''"esc Tdiff 



(4.10) 



the density of CHAMPs in the disc at t > is 



+C exp 



TdiS 
t 



1 + 



Tdc 



exp 



t 

''"dec 



(4.11) 



provided that race ^ T cS (see Appendix || for details). The constant C is fixed by imposing 
the initial condition that n^ d (td) = tiqF'[1 — exp(— t<i/T& cc j\. Since the Galactic disc is 
about 12 Gyr old, we take ~ 1.6 Gyr. 

The number density of positive CHAMPs in the sea water is predicted to be: 



1 

Ad 



"cd 
last 3 Gyr 



K d {t)v cd + nt d {t)v hd )dt 



(4.12) 



where d is the average ocean depth (d ~ 2.6 km), v c d and Vhd are the characteristic velocities 
of particles in the frame corotating with the Sun for the cold and hot population (i.e. v c d — 
17 km s , Vhd ^ 300 km s _1 ). The integration is carried out over the accumulation time 
of the CHAMPs inside the sea water, about the age of oceans, ~ 3 Gyr (e.g., [36, |37|]). 



By requiring that the relative abundance of superheavy isotopes of hydrogen, compared 
to ordinary hydrogen is less than about 6 x 10~ 15 in the range 10 TeV < mx < 6 x 10 
TeV the exclusion diagram was derived for three different values of mx (see Figure 
|2|) for vm = 300 km s _1 . Note that n^ d oc v^ d . For Tdiff we have adopted the maximum 
value permitted by the inequality (O). Those neutraCHAMPs that are converted to X^a 
by charge-exchange scattering when they cross the disc will also contribute as a source of 
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heavy hydrogen in Eq. ( [4,12 ). The inclusion of this additional source of superheavy sea 



water shifts the curves in the top panel of Figure [2| downwards, making the allowed region 
of parameters more restrictive. 

The curves in the T esc versus t^ cc plane at the top panel of Figure |2| define the maximum 
value of t csc compatible with the heavy-water searches, as a function of Td ec , i.e. r™^ x = 
Q(Tdcc, m x)- Along these curves, we have derived the maximum present-day values of n^, 
n^ d and nt d allowed by sea water searches (Figure ||b). We see that in order for CHAMPs 
to be abundant in the halo (say, u^/uq > 0.2), mx > 10 4 TeV is required. It is interesting 
to note that the permitted region of parameters (rd e c>Tesc) is very restricted for a particle 
mass of 2 x 10 3 TeV. 

We see from Figure || that for mx = 2 x 10 4 TeV, T esc < 0.6 Gyr is a guarantee 
that the sea- water constraint is fulfilled. However, according to Eq. Q4.2| ), the population 
of hot CHAMPs should be accelerated to a velocity dispersion of > 3, 400 km s _1 in 
order to escape in less than ~ 1 Gyr. This value is, of course, much larger than our 
reference value v^d = 300 km s _1 used to derive n+ a in Eq. ( 4.12j ). A set of self-consistent 



calculations, which include the fact that v^d and r esc are not independent, shows that 
solutions compatible with the lack of detection of sea water and 300 < v^d < 10, 000 km 
s , require Td ec > 500 Gyr for a particle mass of 2 x 10 4 TeV and r dcc > 2.5 x 10 3 Gyr for 
a mass of 2 x 10 3 TeV. 

It is now clear that the allowed values for Td ec depend on mx but also on the adopted 
value for Vhd- The parameters r esc and Td ec are unconstrained for mx > 6 x 10 4 TeV, 
because the concentration of anomalously heavy hydrogen in sea water is very uncertain. 
In the next section we will consider the vertical pressure equilibrium in the Galactic disc 
and will find out a more robust lower limit on the lifetime of x particles. 

5. The global magnetic support of halo CHAMPs 

In the idealised situation that the magnetic field in the disc, at \z\ < Z m { n , is horizontal and 
the halo is unmagnetized, when charged particles with small gyroradii try to penetrate the 
disc, they execute approximately half a gyro-orbit before finding themselves back in the 
unmagnetized region and with velocities directed away from the magnetized region (Fig. [l]). 
In the boundary layer, a current layer develops as a result of a thermal, unmagnetized 
plasma interacting with a magnetized region. It is a classical result that the (kinetic) 
motions of individual particles in collisionless plasmas can be reconciled with the role 
inferred for the pressure in MHD (e.g., [38, 39]). Since charged CHAMPs in the halo are 



essentially collisionless, the CHAMP momentum flux mx(nt +nT)cr 2 , at the caps of the disc 
should be balanced by the pressure Pb of the magnetic field in the disc. As we discussed 
in the previous section, the galactic disc may contain hot CHAMPs that may participate 
in the pressure balance, whereas particles trapped in the cold phase do not contribute to 
support the halo CHAMPs because its density is very low at \z\ ~ Z m \ n . Ignoring the 
weight of coronal gas and taking = and n^ d = n~^ d , it holds 

PB(Z min ) > 2m x nl<j 2 v - 2m x nl d a 2 hd . (5.1) 
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Figure 2: The allowed regions of r csc versus Td oc parameter space from searches for superheavy 
hydrogen in deep sea water, for different mx (top panel) . At each mass mx , the region of parameters 
below the corresponding curve is allowed. The maximum present-day abundance of neutral and 
ionized superheavy hydrogen, normalized to no, in the halo, cold disc and hot disc, allowed by the 
searches of heavy water in the sea, is shown as a function of Td ec (bottom panel). 



Here is the density at the top of the disc. 

Let us now consider the term n~^ d a\ d in Eq. ( |5.lD Both CHAMPs and cosmic rays in 
the Galactic disc may be subject to Fermi acceleration mechanisms in shock waves. It is 
known that cosmic rays are in approximate energy equipartition with the magnetic field 



in the diffuse interstellar medium [40]. Equipartition arguments are usually adopted to 
find the magnetic field in other external galaxies (e.g., The concept of equipartition 

between the magnetic field and energetic cosmic rays in our Galaxy (or other galaxies) is 
consistent with the widely held belief that the cosmic rays diffuse through the field but do 
not dominate it. In analogy to cosmic rays, we expect equipartition between magnetic fields 
and CHAMPs in the hot disc. If we assume that the pressure exerted by the CHAMPs 
in the hot disc is a fraction of the magnetic pressure, Eq. ([O]) is simplified to Pb > 
2(1 + fihd^n^al, with fihd of order of unity. It is convenient to express this equation in 
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terms of the mass fraction of dark matter in charged CHAMPs in the halo fh as 



P B (Z min ) > fhPovl (5.2) 

where po is the density of dark matter at the top of the disc if it would have not decayed. 
We see that by requiring that the confinement of the magnetic pressure is entirely due to 
the halo charged CHAMPs, an upper value on the abundance of CHAMPs can be derived. 

In a real galaxy, the topology of the magnetic field is more complex. The halo is also 
magnetized and the magnetic field may rise above the disc. In the following, we derive 
a constraint analogous to Eq. (|5.2| ) but for this more general case. Suppose that charged 
halo CHAMPs cannot penetrate down to z = Z mm (with Z m \ n > H) in the lifetime of the 
Galaxy because of the Galactic magnetic barrier (see § [4.1[ ). Integrating the equation of 
vertical equilibrium from z = Z mm to z = oo, assuming that the magnetic field is horizontal 
at Z m [ n and zero pressure at z = oo, we find 

i r°° 

P B (Z mia ) > T— 5- / f h po(R,z)K z dz, (5.3) 
J- + Phd JZ mia 

where po(R, z) is the density of dark matter if it would have not decayed and K z the vertical 
positive gravitational acceleration. Again, Eq. ( |5.3j ) is written as an inequality because the 
weight of coronal gas has been neglected. If the magnetic field is not horizontal at Z m \ n 
because the topology of the magnetic field is such that only in certain areas the magnetic 
field rises above the disc, the effective vertical magnetic pressure can be represented by 
(B 2 - 2B 2 )/8ir with B 2 and B 2 Z interpreted as averages in the (x,y) plane |4~2]j . Therefore, 
the magnetic tension reduces the effective vertical magnetic pressure and the inequality in 
Eq. fl5.3j ) still applies. 

For a spherical dark halo, the weight term of charged halo CHAMPs at the solar 
vicinity is 

oo 

f hPo K z dz = 1.7 f h x l(T 10 dyn cm" 2 

X (o.OlM pc- 3 ) (l50 km s" 1 ) ' ^'^ 

Interestingly, the observed synchrotron emission above the plane in the solar neigh- 
bourhood implies that the scale height of the magnetic field is greater than what would be 
inferred from the weight distribution of the interstellar matter (e.g., |43| ). The observed 
synchrotron emission above the plane in the solar neighbourhood indicates that the total 
magnetic field strength is 2 — 5 p>G at a height of z = 1 kpc [43, 45]. If we identify 



Z m in as the half width at half maximum (HWHM) of the magnetoionic disc ~ 1 kpc (e.g., 
f4q| ) and by evaluating the magnetic pressure at z = Z m - m 1 kpc, we obtain the desired 
constraint on fh, once adopting the highest magnetic value of 5/iG allowed by observations: 



/„< 7X10-(1 + , M )(^L_)- 1 ( IH ^_ T , . (5.5) 
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This estimate is very robust to the precise value adopted for Z m { n because the magnetic 
field decays very slowly with z. 

In our derivation, we have assumed that the halo is spherical. Consider now an oblate 
isothermal dark halo with axis ratio q: 

where v c is the asymptotic circular velocity at the equatorial plane and v = 7 -1 arcsin7, 
with 7 = \J\ — q 2 . In this model, the velocity dispersion is given, within less than 10%, 
by a v ~ 1.16y / g(w c /v / 2), for flattening 0.05 < q < 0.5 (e.g. [47]). Even though the velocity 



dispersion for q < 1 is smaller than in the spherical case, the weight term changes only by 
~ 10% as compared to the spherical case, even for rather flattened haloes (q ~ 0.5). 

Consider now a portion of the disc at larger galactocentric distances, say R = 2R@. 
Following the same procedure as in the solar neighbourhood, we need to estimate the 
total magnetic pressure at (2R & , Z m { n ), which should be responsible to give support to the 
charged halo CHAMPs. The large-scale magnetic field may have a scaleheight 5-10 times 
the scaleheight of the neutral gas disc, so that we may assume that .Bo(^min) — Bq(z = 0). 
The random magnetic field is expected to be roughly in equipartition with the kinetic 
energy in the turbulence. Therefore, its vertical scaleheight should be similar to that of 
the gas. If magnetic fields are still a barrier for halo CHAMPs, then we may assume that 
Z m i Q > H and, consequently, the magnetic pressure by the random component at Z m \ a is 
less than 10% the pressure by the random field at z = 0. Collecting both contributions, 
we derive an upper limit for the total magnetic pressure at Z m [ n : 

P B< M±M! = | (1+ „, a) , (5 . 7) 

where a = b 2 /B 2 , with b 2 and Bq evaluated at z = 0. The ordered magnetic field is 
difficult to measure in the outer Galaxy, but there is evidence that it decays with radius 
R as a power-law between R^ 1 and i? -2 , probably as exp(— R/Rb) with Rb = 8.5 kpc 
I if, @- The uniform magnetic field in the solar neighbourhood is 2-4/iG, depending on 



the authors ]^9|, |50[j . If we generously take a value in the solar circle of 4//G, we infer a 
strength of Bq ~ 1.5^G at 2R & . Assuming a spherical dark halo with a mass density at 
2Rq of ~ po,o/4, then J fhp K z dz = fhPo : Q(T 2 /4. At 2R Q , our assumption that the halo 
is spherical is a very good approximation (e.g., |5l], HH). By imposing pressure balance at 
z = Z m \ D (Eq. |5.3|) , the following constraint for fh is inferred 



fh<2x 10^(1 + 0. la) (1 + /3 M ) 
Po,o/4 



O.OO25M pc- 3 J V 150 km s 



(5.8) 



Other observational estimates assure our generously-taken magnetic intensity. In fact, 
data from rotation measurements of pulsars suggest uniform magnetic fields of ~ 0.7 ^G 
at R = 2Rq f53|]) which coincides with the extrapolation of the fit of radial variation of the 



regular field by Han et al. 4£] 
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Beyond 2R & it is uncertain if supernovae shocks are able to clean the disc from 
CHAMPs. It might be also possible that beyond the optical radius, the magnetic field 
is too weak to prevent CHAMPs from crossing the disc, but any more complicated analysis 
is useless in the face of such ignorance. 

Combining Eqs fl5.5|) and (|5lj|) and taking fihd ~ 1, the present-day fraction of charged 
CHAMPs in the halo must be smaller than (4 — 14) x 1CT 3 , which implies Td cc ^ (0.95 — 
3.4) x 10 3 Gyr. This lower limit on Td ec , which is valid for mj as large as 10 6 TeV, is 
comparable to the constraint inferred from the lack of anomalously heavy water in the sea 
for mx ~ 2 x 10 3 TeV (see section |4~3|) . We conclude that although charged particles can 
be suspended in the halo, so that they would be impossible to detect as they never reach 
the Earth, the mass fraction of charged CHAMPs in the halo must be rather small. 

We must note that, whereas constraints from BBN and heavy- water searches are only 
relevant if CHAMPs are singly charged, because for other charges, a CHAMP no longer 
behaves as a proton, the constraint fh < (0.4 — 1.4) x 10~ 2 derived in this section, is 
independent of charge e, as long as charged CHAMPs and magnetic fields are in pressure 
equipartition in the Galactic disc. 



6. Ram pressure stripping and collisions of galaxy clusters 

Magnetic fields couple charged CHAMPs with themselves and with ordinary matter. This 
coupling might cause ram pressure stripping of both baryonic and DM of subhaloes and 
satellite systems. Consider, for instance, the collision of two galaxy clusters. Estimates 
for the magnetic field strength in clusters range from roughly 1 — 10 fiG at the center and 
0.1 — 1 /iG at a radius of 1 Mpc. With these values, the ratio between thermal pressure 
Pth and magnetic pressure for the CHAMPs is j3 = SirPth/B 2 m 2/ x 10 3-4 , that is, a hot 
plasma. Even in this dynamically weak magnetic field, the mean gyroradius for a CHAMP 
with mx = 10 6 e TeV, is < 5 pc at the center and < 50 pc at 1 Mpc. The governing 



equations of collisionless hot plasmas were developed by Chew et al. [54], whose theory is 
known as the Chew-Goldberger-Low approximation. This approximation, which leads to 
MHD equations with anisotropic pressure, is satisfactory when the Larmor frequency is 
large compared to other characteristic frequencies of the problem and the mean particle 
gyroradius is short compared to the distance over which all the macroscopic quantities 
change appreciably (e.g., |)5|, |5q| ). Therefore, charged massive particles in the halo of 
galaxy clusters can be described in the fluid-like anisotropic MHD approximation; in the 
merger process, they would behave as a clump of fluid, experiencing ram pressure stripping 
and drag deceleration similar to the gas component. Since CHAMPs should be attached to 
the gas component, the observed offset between the centroid of DM and the collisional gas 



of the subcluster in the Bullet Cluster implies / < 1 (e.g., (57), |58|). Although the current 
lensing data accuracy is not sufficient to derive the mass distribution of the subcluster in 
the Bullet Cluster, the derived mass estimates of the subcluster leave little room for DM 
in the gas bullet. 

Galactic halo CHAMPs may also exert ram pressure on the gas component of the LMC 
and its stream due to their continuous scattering by the intrinsic magnetic field of the LMC 
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and the Magellanic stream. For a Milky Way- type halo of ~ 10 12 M Q , a fraction of 
4 x 10~ 3 implies that the mass in charged CHAMPs could be up to ~ 4 x 10 9 M and the 
density at 50 kpc of 1.2 x 10~ 6 M pc~ 3 . Since these values are smaller than those required 
to explain the mass and extension of the Magellanic Stream and the size and morphology 
of the gaseous disc of LMC [59|, we cannot reduce any further our upper limit on fh with 
the current observations of the LMC disc and the Magellanic Stream. 



7. Conclusions 

Whilst the common wisdom holds that DM is neutral and collisionless, it is important to 
explore the possibility of it having nonzero, not necessarily integer, charge. If a fraction of 
the mass of haloes is made up by charged CHAMPs, it may have a strong impact on the 
observable Universe because of the coupling between magnetic fields and CHAMPs. For 
instance, ejection of charged CHAMPs from the regions with intense magnetic fields, i.e., 
from the central parts of galaxies, would help alleviate the cuspy halo problem. In this 
work, we have constrained the present abundance of CHAMPs in galactic haloes. 

We have explored a model where neutral dark matter decays into non-relativistic 
charged products. From BBN and CMB, we find that the decay lifetime should be > 0.1 
Gyr. The non-detection of heavy sea-water puts a limit on the timescale for charged 
CHAMPs to escape from the Galactic disc. We have considered the pressure support of 
CHAMPs in our Galaxy to derive a simple, upper limit on the fraction of CHAMPs and 
milliCHAMPs in galactic haloes. Assuming that the accelerated CHAMPs in the disc are 
in pressure equilibrium with the magnetic field, we find that fh < (0.4 — 1.4) x 10~ 2 . This 
constraint rules out CHAMPs as the origin of the cores in LSB and dwarf galaxies. The 
reduction of the central density after they have driven the formation of galactic haloes 
would be insignificant. Even if all the CHAMPs were depleted from the central parts of 
the galaxies, the rotation velocity in a certain galaxy would suffer a negligible change of 
(0.2 - 0.7)% for f h ~ (0.4 - 1.4) x 10~ 2 . In the range of astrophysical interest, CHAMPs 
behave like strongly interacting (fluid-like) dark matter (SIDM). Thus, they face many of 
the problems attributed to SIDM. As some examples, we have discussed the survival of the 
Magellanic Stream and the mass distribution of the Bullet Cluster. Our constraint that 
the mass in CHAMPs in the Galaxy is not larger than the mass of coronal gas in the halo 
seems to apply also to galaxy clusters. 
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A. Recombination and ionization of CHAMPs in the Galactic corona 



Positive CHAMPs can recombine with free electrons in the Galactic corona to form neutral 
bound atoms (X^e~). As X±s behave exactly like protons with thermal velocity disper- 
sion 4 of 150 km s _1 , the fraction of nuclei that are neutral is comparable, or even 
smaller, than the neutral fraction of hydrogen in the coronal gas, which is extremelly small 
(~ 10~ 6 , e.g., and therefore it can be ignored. 

The recombination rate of negative CHAMPs with protons in the Galactic halo includ- 
ing the ground n = 1 level is given by a^n p n^, where the coefficient for recombination 
is: 

r 71=1 

where v 2 s = v\ + v 2 (e.g., |j, At the Galactic halo, we have v e g ~ 180 km s _1 and 

thus 0$ = 4 x 10" 18 cm 3 s _1 . 

Ionization of (X^p) by collisions with and Xf are also important for the determi- 
nation of the abundance of neutraCHAMPs in galactic haloes. The coefficient for collisional 
ionization is: 



7i< 



1.3 x W-'T^FE^eY) exp [cm 3 /s], (A.2) 



where -Ebin — 25 keV is the binding energy of the atom (X^p), Tx the temperature of bare 
CHAMPs and F ~ 0.83 for hydrogenic atoms (e.g., f62|). For massive CHAMPs in the 
Galactic halo, we have 

-io ( m x \ 1/2 ( 



7ion = 1.0 x 10- 1U — A 7 ; [cm 3 /s]. (A.3) 

\2xl0 4 TeVy V 150km / S / 



Putting together, 



dn. 



+ lion(n h +nj)n nC - a y R 'n p n h , (A.4) 
at r dcc 

where n n c = — is the density of neutraCHAMPs in the halo. The abundance of 
neutraCHAMPs as compared to CHAMPs at ionization equilibrium is ~ 1.0 x 10~ 3 . This 
estimate should be considered as an upper limit because (X^p) may be converted to a 
charged state, (X^a), by a charge exchange reaction. Unfortunalely, the exchange cross 
section is very uncertain Q. 

Coronal neutraCHAMPs can penetrate the Galactic disc, reach Earth and stop in the 
atmosphere or ocean Searches for coronal neutraCHAMPs in cosmic rays rule out 
particles with masses between 100 and a few 10 4 TeV if all the X~ are bound to a proton 
and the charge exchange cross section with C and O nuclei is in the interval from 30 mb 
to 30 b p|,|6l M. 



4 



CHAMPs and protons have different temperatures. The thermal (one-dimensional) velocity dispersion 
of protons of the hot gas (~ 10 6 K) at the Galactic corona is v p = yJk B T/m p « 100 km s" 1 , which is a bit 
smaller than the adopted value for the velocity dispersion of dark matter particles in the halo (~ 150 km 
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B. Solving the differential equation for n 



hd 



In section 4.3, we derived the differential equation for the number density of hot X + in the 
disc as: 



dn 



hd 



dt 



n hd , n h n hd 



ec ' esc 



Tdiff 



with 



and 



n x (t) = n exp ( —J , 

V r dec/ 



n h (*) = n o 
Inserting them into Eq. (|B.1| ), we find 



1 — exp 



t 

''"dec 



dt r eS r diff 
where the definition of r e g was given in Eq. (pLlOD and 



5 (t) = 1 + F' 



1 exp 



The general solution of Eq. 



> T diff 
' dec 



when Tdec ^ r eff is 



(-—) f/^(*)exp f-Vt + C 

V r eff/ U r diff \ r eff/ 



n hd(*) = exp - — 



where C is a constant. The integral can be performed analytically: 



t 



Ted 



exp / g{t) exp — = r eff 



1 + 



^Tdiff - 7dc 
Tdec — Teff 



exp 



(B.l) 
(B.2) 
(B.3) 

(B.4) 
(B.5) 
(B.6) 
(B.7) 



The resulting form for n~^ d (t) is given in Eq. ( [4.11 ). 
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